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Nomenclature 


d = diameter of the suction hole 

m = - Ps^s/^eUe, suction ratio 

Re 0 = Reynolds number based on the momentum 

thickness and the edge velocity 

U,W,V = mean velocity components along x-, y- and 

z-directions, respectively 

U e ,W e = mean velocity components at the edge of the 

boundary layer 

Vi,V 2 ,V 3 = velocity components normal to the axis of beam- 

pairs 1, 2, and 3, respectively 

V s = suction velocity 

u, v , w = fluctuating velocity components along x-, y- and 

z-directions, respectively 

u 2 } v 2 w 2 = Reynolds normal-stress components along x-, y- 

and z-directions, respectively 

Reynolds shear-stress components along x-, y- 
and z-directions, respectively 

friction velocity 

streamwise, normal and transverse coordinates 
flow turning angle 

cross-flow angle between the external and local 
streamlines 

P w = cross-flow angle, measured as the angle between 

the projected external streamline and the surface 
shear stress vector (limiting streamline) 
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5 


= boundary layer thickness 


4> = inclination of the laser-optics table to the 

horizontal plane 

v = kinematic viscosity 

P Sj P e = density of suction- and freestream air, 

respectively 

0 = inclination of the beam-pair axis to the flow axis 

Q x , Qy, Q z = mean vorticity components along x-, y- and z- 

directions, respectively 
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Abstract 


This report provides an over view of the three-dimensional 
turbulent boundary layer concepts and, of the currently available 
experimental information for their turbulence modeling. It is found 
that more reliable turbulence data, especially of the Reynolds stress 
transport terms, is needed to improve the existing modeling 
capabilities. An experiment is proposed to study the three- 
dimensional boundary layer formed by a "sink flow" in a fully 
developed two-dimensional turbulent boundary layer. Also, the 
mean and turbulence field measurement procedure using a three- 
component laser Doppler velocimeter is described. 

1. Introduction 

By definition, a three-dimensional boundary layer has three 
non-zero mean velocity components that are functions of the three 

reference coordinates (x, y and z). The external potential velocity is 
a function of two coordinates (x, z) in the plane of the surface. If the 
external potential streamlines are straight lines which either 
converge or diverge then, as compared to a two-dimensional 
boundary layer there is only change in the boundary layer thickness 
and the velocity vector points in one direction throughout the 
boundary layer. On the other hand, if the potential streamlines are 
curved then, the imbalance between the decreasing centrifugal force 
(due to velocity reduction in the boundary layer) and the impressed 
radial pressure gradient force gives rise to secondary flow in the 

boundary layer with skewed velocity profiles as shown in Figure 1. 

This type of secondary flow is called by various names, such as, skew 
induced secondary flow (Prandtl’s first kind of secondary flow- 
Bradshaw 1987), pressure-driven secondary flow (Anderson and 
Eaton 1989, Baskaran et al. 1990 etc.). In Figure 1, the x axis is 
chosen along the direction of the potential streamlines, z axis is 
orthogonal to x axis in the surface plane and y axis is chosen normal 
to the surface. Accordingly, the primary (U) component of velocity is 
the x-component, vertical (V) velocity is the y-component and 
secondary or cross-flow (W) component of the velocity is the z- 
component. Within the boundary layer centrifugal, viscous and 

convection effects are significant, but in the potential region only the 
centrifugal and convection effects are important. Therefore, if the 
sign of the lateral pressure gradient changes along external 
streamlines the corresponding cross-flow velocity changes tends to 
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lag behind and gives rise to a S-shaped cross-flow profiles some 
distance downstream of pressure gradient sign change. The velocity 
vectors within the boundary layer are non-collateral and therefore, 

the streamlines are skewed towards the concave side of the potential 

streamlines. The streamlines on the surface (also known as the 

limiting streamlines) have the maximum skew angles and the skew 
angle, in general, is a function of the y-distance normal to the 

surface. Also, the streamwise (3U/3y) and the cross-flow (3W/3y } 
velocity gradients give rise to mean streamwise vortices. Once this 
type of skew induced vortices are formed they are further diffused 
by viscous and Reynolds stresses. 

There is a second kind of three-dimensional boundary layer 
that is typified by a turbulent flow in the rectangular corner or of a 
wing-body junction. Here, the streamwise mean vortices are created 
by shear stress gradients in both the y- and z- directions. These 
type of secondary flows are referred to as "stress induced secondary 
flow" or "secondary flow of Prandtl's second kind". Obviously, these 
vortices occur only in turbulent flows and, are, generally, weaker 
compared to the skew induced vortices. 

2. Three-dimensional Boundary Layer Concepts 

In three-dimensions generally we use a coordinate system that 
conforms to the surface (Body Fitted Coordinates), and this involves a 
network of two families of lines, not necessarily orthogonal, mapped 
over the surface plus the normals to the surface. Hence the 
curvature of the coordinates are introduced with associated 
centrifugal and coriolis terms in the equation as well as metrics of 
the system. 

However, for small surface curvature, cartesian coordinates can 
be used for convenience. The boundary layer equations for the 
incompressible flow are given by 

_ DU 3p 

D = + 

Dt 3x 3y 

DW 3p 3x 

K Dt 3z dy 

au av aw 

3x 3y 3z 
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where, x xy = p 


au 

dy 


, aw 

an dT„ = H — 


for laminar flows 


and, x = p— — puv and x = ]LL— pwv for turbulent flows 

dy dy 

dp/dy is assumed negligible and, dp/dx and dp/dz are known from the 
potential flow solution. In turbulent flows, the U-component motion 
is governed by d uv /dy and the W- component is governed by 

a vw /ay. in certain flows the term d^ w /dx, omitted from the above 
equation, may also be significant. The three-dimensional Navier- 
Stokes equation are elliptic, which means that information from a 
point can be propagated by convection, viscous diffusion, or pressure 
variations. Under the three-dimensional boundary layer 
approximations the influence of the solution at one point is 

transfered to other downstream points within a wedge-like region 
formed by the surface streamline (actually surface shear stress 
vector) and the free streamline by viscous diffusion in the vertical 

direction and by convection in the streamwise/cross-stream 
directions. This region is called "domain of influence" as indicated in 

the Figure 2. On the other hand, the solution at a given point 

depends on the solution in another upstream wedge-like region 
called the "domain of dependence". Also, in three-dimensional 
boundary layers streamwise vortices, Q x are generated, as given by 
the following vorticity transport equation, by; (i) quasi-inviscid 
deflection of existing mean vorticity Q y and Q z by the mean velocity 
gradients and, (ii) turbulent stress gradients. 


U 




dx 


V^l + W 
3y 


dc\ _ au _ au _ au 

a 2 rr — t 


f 


i2 X 


V 


dy 2 dz 2 


dydz 


(v 2 - w 2 ) + vV 2 ^ 


(-vw) + 

( 1 ) 


For the purposes of analysis, three-dimensional boundary layers are 
classified in to three categories; thin shear layers, slender shear flow , 
and three-dimensional shear flow. 

2.1 Thin shear layer 

In thin shear boundary layer flows, d/dy » d/dx ~ d/dz when 
operating on any velocity component. These correspond closely to 
two-dimensional boundary layers and sometimes called boundary 
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sheets. Here, V is small compared to U and W and the pressure is 
constant in the y-direction. Such flows occur on yawed flat plate, 
finite wings and propeller blades, except in the region near the root 
and tip. Three-dimensional thin shear layers have been solved in 
specific cases using some further simplifying assumptions. These 
are; (i) the "independence principle" according to Prandtl (1946), (ii) 
the infinite swept wing or infinite yawed cylinder and (iii) the "small 
cross flow assumptions" or the "principle of prevalence". 

In the case of thin shear layers, with negligible viscous 
diffusion in the x- and z-directions the information is propagated in 
the cross-flow plane only by convection at an angle tan-!(W/U). 
According to the independence principle the velocity in the boundary 
layer which is parallel to the wall is also parallel the potential flow at 
all points, or in other words, W/U= W e /U e . This is valid for a yawed 
flat plate boundary layer. This causes the boundary layer thickness 
of a turbulent boundary layer on a yawed flat plate (e.g. Ashkenas 
and Riddel 1955) grow faster in the downstream direction than with 
an unyawed plate. But, the independence principle can be applied to 
laminar flows only, because in turbulent flows there is a strong 
interaction between the spanwise and chordwise components of the 
velocity fluctuations which change the mean velocities. 

For infinite swept wing or infinite yawed cylinder conditions, 
the z-derivative is zero and the surface is assumed to be developable 
(one that can be obtained by rolling up a plane flexible sheet). The 
calculation methods for this case is provided by Adams (1975), 
Krause (1974) and Cebeci (1974). 

The "small cross-flow" or the "principle of prevalence" 
assumption amounts to assuming that W is so small that it may be 
ignored in the momentum equation for U. Generally, small cross- 
flow assumption holds good in turbulent boundary layers to cross- 
flow angles, (3 W smaller than six degrees. In flows which have a 
plane of symmetry some further simplifications are possible with 
small cross-flow assumption. For example, Nash and Patel(1972) 
assumed the cross-flow mean velocity as well as the shear stress, 

pvw, to vanish whereas the gradients in the z-direction are retained 
for their three-dimensional boundary layer with a plane of 
symmetry. A more detailed discussion of the turbulent flow at a 
plane of symmetry are available in Johnston (1960) and 
Pierce(1963). 

2.2 Slender shear flow 

This is the kind of shear flow along corner walls, where 8/3 y, 
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d/d z » d/d x. The shear stress gradients in both y- and z-directions 

are important in the Reynolds stress equations. This can give rise to 
discrete streamwise vortices in the boundary layer as described 
above (c/. equation 1). 

2.3 Three-dimensional shear flow 

The third category is called three-dimensional shear flow, 
where d/dx ~ d/dy ~ d/dz. This needs higher order boundary layer 
equations with term-by-term modeling of Reynolds stresses, or the 
full Navier-Stokes equations. 

Detailed reviews of all the above three kinds of three- 
dimensional flows may be found in; Cooke and Hall (1962), Joubert et 
al. (1967), Sherman (1968), Wheeler and Johnston (1972,1973), Nash 
and Patel (1972), Blottner (1975), Horlock et al. (1966), Fernholz 
(1982) and Eichelbrenner and Oudart (1955). 

3. Previous experiments in three-dimensional turbulent 

boundary layer flows 

A few experimental results are available for three-dimensional 
thin shear layers: on pressure-driven secondary flows on swept 
wings by Etheridge (1971), van den Berg et al. (1975) and Elsenaar 
and Boelsma (1974); on curved channel flows by Klinksiek and Pierce 
(1970), Vermeulen (1971) and ship hulls by Larsson (1975). 
Experiments on shear-driven secondary flows are even fewer; 
Bradshaw and Terrell (1969), Crabbe (1971) and Driver and Hebbar 
(1985). But in all the measurements, no accurate turbulence data in 
the vicinity of the wall is available. 

3.1 Mean flow data 

Bradshaw (1987) considered the effects of imposing a spanwise 
pressure gradient on a two-dimensional boundary layer that lead to 
a three-dimensional flow. Here, the spanwise pressure gradient 

(albeit small) is created as a result of changes in the radius of 
curvature of streamlines as the flow passed over a curved section. 
The streamlines at radii greater than the mean radius experience 
lower streamwise velocities and the streamlines at radii smaller than 
the mean radius experience higher streamwise velocities. This 
spanwise velocity gradient (dU/dz) gives rise to a streamwise 
vorticity component, Q x , by laterally skewing the pre-existing 
spanwise vorticity vector, D z (predominantly -dV/dy). The axial 
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vorticity created this way has mostly contributions from 3W/3y as 
dV/dz is significantly small, at least initially. For small flow 
deflection angles a, we have the Squire Winter Hawthorn (SWH) 
secondary flow formula, d(£2 x /£2 z )/dx = d(W/U)/dx or 2a = - tan' 
HW/U)), which implies that vortex lines are skewed through an 
angle equal and opposite to that through which the flow has turned. 
This relation is found to agree reasonably well with the experiments 
(e.g. Bradshaw and Pontikos 1985, Driver and Hebbar 1985) for the 
external streamlines. However, close to the surface the flow 

deflection angles are smaller than SWH predictions. According to 
Bradshaw (1987), when the sign of the spanwise pressure gradient 
changes further downstream a crossover profile in W component 
may occur. It is difficult to fit crossover profiles by means of simple 
relations. For non-crossover profiles Johnston (1960) proposed a 
correlation model for a composite velocity defined as U/cos(3 w as, 


U 

u * COS p w 



+ C 


Here, u* is the friction velocity formed with total shear stress, and P w 
is the difference between the direction of the velocity vector in at 
the outer edge of the boundary layer and the skin friction vector at 
the wall. Among several other mean velocity models for Mager's 
(1952) relation, 



seems to provide the best agreement with experiments. 

3.2 Turbulence data 

Turbulence data are available for a variety of three- 
dimensional boundary layer flows including pressure-driven 
secondary flow, vortex imbedded boundary layer flow and cross- 
stream wall jets . These results are discussed now. The turbulence 
results of Bradshaw and Pontikos (1985), from their study on infinite 
swept wing, indicate some interesting behavior: (i) the ratio of the 
magnitude of the resultant turbulent shear stress to the turbulent 
kinetic energy drops rapidly as the cross flow velocity component 

increases, (ii) the development of the secondary shear stress, pvw is 
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slower than the development of cross-flow velocity W and its 
gradient 8W/8y, and (iii) the viscous dissipation is larger than would 
be for a two-dimensional flow. Overall, the turbulent activity was 
decreased by the cross-flow, leading to reductions in turbulent 
transport of momentum, turbulent energy and turbulent shear stress 
across the layer. The explanation is that the large eddies in the 

initial two-dimensional flow are tilted sideways by the cross-flow 
velocity gradient, and this tilting reduces the capacity of the eddies 
to transport the above three turbulent quantities in the transverse 
direction. However, the pressure-strain" redistribution terms and 
the turbulent transport terms in the Reynolds-stress transport 
equations may still respond immediately to the 3U/9y and 9W/9y 
gradients. Anderson and Eaton (1986) also noticed similar effects on 
flow around surface mounted obstacles 

In the 1982 Eurovisc Workshop, existing turbulent models and 
calculation methods have been evaluated for predicting flow on 
infinite swept wing. The general consensus is that no single model or 
calculation methods provide satisfactory predictions of all the flow 
quantities, particularly in the near wall regions. It is also deliberated 
that detailed turbulence measurements in the sublayer are scarce 
and models have to be adjusted to yield the right mean velocity 
profile. 

If the lateral deflection of the vortex line is confined to only 
small spanwise distance, then concentrated longitudinal vortices are 
formed. For example, the trailing vortices created by the leading 
edge of a slender wing. These vortices can be isolated vortices (e.g. 
Gad-el-Hak and Blackwelder 1985, Cutler and Bradshaw 1986) in 
inviscid flow or they can be embedded in the boundary layer (e. g. 
Shabaka et al. 1985, Westphal et al. 1985, Eibeck and Eaton 1985, 
Mehta and Bradshaw 1986 and Subramanian et al. 1993). Cutler and 
Bradshaw's (1986) studies show that in the core of the isolated 
vortex the turbulent mixing is small yet, longitudinal turbulence 
intensity values are significant due to the presence of longitudinal 
waves. Since turbulence production per say is absent, a Reynolds 
stress model which works well for truly turbulent flows may not 
work for this case (e.g. Majumdar and Rodi 1985). 

Concentrated vortices from an oncoming boundary layer are 
generated in junction flows of a wing, a turbomachinery blade, or a 
building. If the leading edge is sharp, the concentrated vortices may 
be close enough to surface to be rapidly diffused by viscous or 
Reynolds stresses. A detailed investigation of turbine-blade junction 
flow has been done by Langston et al. (1977) and Langston (1980), 


but complete turbulence data is not available. Shabaka and 
Bradshaw (1981) have reported results of turbulence measurements 
on a simplified wing body junction flow. Here one can evidence the 
formation of pairs of stress induced vortices following the initial 
formation of pressure-driven horse-shoe vortices, as shown by 
Nakayama and Rahai (1984). 

The general consensus is that for both corner vortex as well as 
embedded vortex flows even the most refined turbulence models do 
not give adequate predictions of the cross-stream intensities and 
secondary shear stresses that control the diffusion of streamwise 
vorticity. Probably, an improved modeling of the pressure-strain 
term in the Reynolds transport equation is necessary. 

Flows in curved channels and ducts develop secondary flow 
(and longitudinal vortices) due to centrifugal instability. Although 
these vortices are weaker than the junction- and embedded vortices 
(Subramanian et al. 1992) they can change the turbulence structure. 
Longitudinal surface curvature tends to stabilize the turbulence near 
the inner wall and destabilize that near the outer wall. The most 
recent detailed turbulent measurements are by Chang et al. (1983) in 
a squared curved duct and Azzola and Humphrey (1984) in a curved 
pipe. But these measurements do no include secondary shear stress 
or triple product data. Unfortunately, the yz-plane shear stress, pvw ? 
is an essential part of the streamwise vorticity generation ( cf . 
equation 1). 

The most spectacular stress induced secondary flows are found 
in three-dimensional jets and wall jets. In this case there is a strong 
interaction between the mean flow and the turbulent stresses. As 
shown by Launder and Rodi (1983) for a three-dimensional wall jet, 
stress induced vortices play as much important role as the skew- 
induced vortices in the dynamics of motion. 

4. Proposed experiment: On the structure of turbulence in a "sink 

flow" 

The above review suggests that for modeling three- 
dimensional turbulent boundary layers flows of engineering interest, 
we still lack the complete understanding of the relation between the 
three-dimensional mean flow and the turbulence structure. The 

main question is the behavior of vw due to the mean flow 
interaction. Conventional measurements of vw (e.g. Andreopolous 
and Rodi 1984, Shayesteh et al. 1985) using hot-wire anemometry 
have calibration problems due to subtraction of two cross-wire 
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readings. Reynolds transport equation for this shear stress 
component involve pressure strain-rate product term which is also 
difficult to measure with reliable accuracy. The proposed 
experiment is aimed at a thorough investigation of mean and 
turbulent fields and their interactions in a boundary layer flow using 
a three-component laser Doppler velocimeter. 

4,1 Experimental setup 

The experiments are to be conducted in the 20" X 28" Shear 
Flow Wind Tunnel facility of Experimental Flow Physics branch at the 
NASA Langley Research Center (LaRC). This tunnel has a rectangular 
test section, 0.51 m X 0.71 m, and 4.57 m long and can produce 
velocities of up to 46 m/s. At a free stream speed of 20 m/s the 

freestream turbulence intensity is less than 0.08 percent 

(Subramanian 1992). A splitter plate with a semi-elliptic leading 
edge can be used as the test surface. The tunnel pressure gradient 
can be adjusted by adjusting the tunnel upper wall. 

A three-dimensional boundary layer can be created by 
introducing suction at a streamwise location x = 3.67 m from the 

leading edge in a nominally two-dimensional turbulent flat plate 
boundary layer. The experiments are to be performed at a Reynolds 
number based on the momentum thickness, Ree, of about 8100 so 
that the turbulence is fully developed and has broad band spectrum 

of distinguishable length scales. Boundary layer suction can created 
by drawing air with a suction/vacuum pump through a 1 cm 
diameter hole (d) at the center line of the plate. At the location of 
the hole the canonical boundary layer thickness, 8, is estimated to be 
about 6.2 cm. Thus the ratio d/8 will be about 0.16. The magnitude 
of the suction can be varied so that the ratio of suction mass flux to 
the free stream mass flux, - Ps^s/Pe^e, called suction ratio m, is in 
the range 0.5-1. 5. For this range the flow through the suction duct 
should be laminar for U e = 20 m/s. At lower suction ratios, only the 
near-wall region of the boundary layer will probably be affected, but 
at higher ratios the entire boundary layer is expected to be altered. 
The "sink flow" created by the suction is meant to introduce the 
additional, V and W velocity components in the flow, particularly 
near the wall. The three mean velocity components, six Reynolds 
stresses and ten Reynolds triple products (contributing to turbulent 
transport) can be measured using a three-component laser Doppler 
velocimeter described below. These measurements should provide 
us an insight as to how the turbulence structure (especially in the 
near-wall region) changes in a sink flow. 
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4.2 LDV measurement 

The LDV system : As mentioned above, for measuring velocities a 

three-component, Argon-ion laser Doppler velocimeter is to be used. 
A LDV system has been developed and being used by Dr. J. T. 
Kegelman and his associates at EFPB. This system consists of a 
COHERENT Innova 90, 5-Watt Argon-ion, water-cooled laser, THERMO 
SYSTEMS, INC., (TSI) transmitting and receiving (photo-multiplier 
tubes) optics for violet (476.5 nm), blue-green (496.5 nm) and green 
(514.5 nm) colored beams, TSI frequency shifter and MACRODYNE, 
INC. signal processing unit. For traversing the probe volume, a 
KLINGER three-axis micro-traverse control that allows micrometer 
movement, is used. The signal processing and traversing can be 
automated using the AEROMETRICS Doppler Signal Analyzer 
software. 

In any turbulence measurement the probe resolution and 
accuracy are important factors. By using a 600 mm lens with a beam 
separation of 50 mm we can obtain an average probe-volume 
diameter of about 120 pm and a length of about 800 pm. Although, a 
subminiature hot wire probe perhaps can provide a better spatial 
resolution, the LDV is preferable from the point of view of non- 
intrusive probing capabilities and not requiring a precalibration. But 
several past studies (e. g. Driver and Hebbar 1985) show that LDV 
measurements are not devoid of problems. For example, non- 
uniform particle seeding can not only slow down the data rate but 
also can give rise to other problems such as multiple seeds in the 
probe-volume at a given time and signal cross-talks (for 2- and 3- 
component LDV). Also, if the beam-pair crossings are non-coincident 
at the probe-volume that will cause velocity bias errors. The errors 
due to these problems can be minimized by careful beam alignment 
and by controlling particle seeding. For air, 2 pm diameter 
polystyrene latex microspheres are good for seeding. In the setup 
shown in Figure 4, the three different color beam pairs intersect 
(non-orthogonally) at the same point in space. The beam pair 1 
(voilet) is oriented at - 55 (= 6 ]) degrees to the flow axis, the beam pair 
2 (blue-green) normal to the flow axis and the beam pair 3 (green) at 
+ 55 (= 03 ) degrees to the flow axis. Then, simultaneous measurements 
of all three velocities can be made. The measured velocities (Vj, V 2 , 
V 3 ) with respect to the beam axes are then transformed into 
cartesian velocity components (U, V, W) before performing the 
ensemble averaging. The transformation matrix is, 



"V 


"u" 

cos0,cos<|> 

0 

sinOjCOSt)) 

v 2 

= 

V 

0 

COS(J) 

0 

V,_ 


w 

COS0 3 COS<t) 

0 

sin0 3 cos(|) 


The optics table has a small dip angle of about 10 (=<{>) degrees with 
respect to horizontal plane in order to facilitate measurements very 
close to the wall. 

Velocity Survey Map: Velocity measurements are proposed around 

the suction hole over an area covering 1.5 diameters upstream to 1.5 
diameters downstream from the hole center and, 1.5 diameters on 
either sides from the hole center (see Figure 5). An uniform traverse 
step size of 0.2 diameter can be used in both x- and z- directions, 
amounting to 225 probe locations per horizontal plane. In the 
normal direction, about 40 non-uniformly spaced heights can be used 
to cover the entire boundary layer. The smallest height 
corresponding to the probe-volume semi-diameter is about three 
viscous units at a freestream velocity of 20 m/s. Therefore, we can 
probe even within the linear sublayer. 

5. Summary 

This report provides an over view of the three-dimensional 
turbulent boundary layer concepts and the existing turbulent 
modeling information. It is concluded that the existing turbulence 
models fall short of accurately predicting boundary layers with 
strong secondary flows, especially when turbulent stress interactions 
also occur. In particular, to model the generation of shear stress 
induced streamwise vortices we need a better understanding of the 
behavior of Reynolds stresses and pressure strain-rates in three- 
dimensional flows. Here, a method is suggested using suction to 
create a three-dimensional boundary layer in a fully developed two- 
dimensional turbulent boundary layer. Since the surface boundary 
condition is altered, the changes to the turbulence structure 
corresponding to the changes in the mean motion is expected to be 
immediate. A non-intrusive, mean and turbulence field 
measurement using a three-component laser Doppler anemometer is 
prescribed. 
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Figure 1. Typical velocity profile shape in a pressure-driven three 

dimensional boundary layer. 
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Figure 2. Regions of influence and dependence in a three- 
dimensional boundary layer. 
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Figure 3. Gruschwitz-Johnston polar plot of W vs. U in streamline 

coordinates. 
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Figure 4. 



Figure 5. Velocity survey mapping region. 
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